Andreev-Bashkin effect and knot solitons in interacting mixture of a charged and a 
neutral superfuids with possible releveance for neutron stars 
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We discuss a mixture of interacting a neutral and a charged Bose condensates, which is supposed 
being realized in interior of neutron stars in form of coexistent neutron superfluid and protonic 
superconductor. We show that in this system, besides ordinary vortices of S^ map, the 

neutron condensate also allows for (meta)stable finite-length knotted solitons, which are 

characterized by a nontrivial Hopf invariant and in some circumsatnces are stablized by Faddeev- 
Skyrme term induced by drag effect. We also consider a helical protonic fluxtube in this system and 
show that, in contrast, it does not induce a Faddeev-Skyrme term. 



INTRODUCTION 



In a standard model for a neutron star its interior fea- 
tures superfluidity of neutron Cooper pairs and supercon- 
ductivity of proton Cooper pairs (see e.g. Both 
these condensates allow vortices of — + map. Ear- 
lier it was suggested that the phenomenon of glitches in 
Crab and Vela pulsars is connected with vortex matter in 
these stars 0. This remains a topic of intensive studies 
and discussions (for recent developments and citations 
see 0). Besides that a standard model for a neutron 
star is a special system being a mixture of interacting a 
charged and a neutral condensates which makes it also 
being a topic of abstract academic interest since such 
a system allows for interesting phenomena with no di- 
rect counterparts in e.g. superconducting metals. Stud- 
ies of topological defects in a mixture of a charged and 
a neutral condensates, so far, concerned only ordinary 
Abrikosov-like columnar vortices (see e.g. |^ Q and ref- 
erences therein). In this paper we argue that, possibly, 
this is not the only one type of stable topological defects 
allowed in neutron stars. We show that due to the drag 
effect in a mixture of a neutral and a charged superfluid 
(Andreev-Bashkin effect) the system also allows under 
certain conditions stable finite-length topological defects 
characterized by a nontrivial Hopf invariant, more pre- 
cisely a special version of knot solitons. 

Finite-length topological defects characterized by a 
nontrivial Hopf invariant were attracting interest for a 
long time in condensed matter physics: earlier it was dis- 
cussed in spin-1 neutral superfluids [1,13, in magnets [lo| . 
in charged and neutral two-component Bose condensates 
[l3| . in spin-triplet superconductors Q and 
in other systems. In neutral systems finite-length closed 
vortices are not stable against shrinkage unless their size 
is stabilized by a conservation of some dynamic quantity. 
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like in case of a propagating vortex loop. A special case 
is a neutral two-component system with a phase separa- 
tion, where a vortex loop made up of one condensate with 
confined in its core circulating second condensate is sta- 
ble against shrinkage • Intrinsically stable topological 
defects characterized by a nontrivial Hopf invariant (the 
knot solitons) have been discussed in the Faddeev nonlin- 
ear 0(3) sigma model |0 where its stability is ensured 
by a special fourth-order derivative term: 

Ff = (diif + a{n-dinxdjnf + K{l-n-na)^, (1) 

where n — (ni, n2,n^) is a three-component unit vector. 
A knot soliton (being in the simplest case a toroidal vor- 
tex loop) is a configuration where the vector n resides 
in the core on e.g. the south pole of the unit sphere, at 
infinity it reaches the north pole, while in between the 
core and the vortex boundary it performs n rotations if 
one goes once around the core and m rotations if one 
goes once along a closed curve in toroidal direction. The 
stability of knots in this model was extensively studied in 
numerical simulations 0| . Recently it was realized that 
this model is relevant for wide class of physical systems. 
First, it was suggested that this model may be relevant in 
the infrared limit of QCD with the knots solitons being 
a candidate for glueballs Besides that an extended 
version of Faddeev model has been derived for two-band 
superconductors ^-nd for triplet superconductors 0| . 

Below we discuss a possibility of formation of finite 
length stable topological defects in a mixture of interact- 
ing charged and neutral Bose condensates. 



A MIXTURE OF INTERACTING 
CONDENSATES 



A mixture of a charged (made up of protonic Cooper 
pairs) and neutral (made up of neutronic Cooper pairs) 
Bose condensates in the interior of neutron stars can 
be described in the hydrodynamic limit by the follow- 



2 



ing Ginzburg-Landau functional 



F = + \r'< + r^v ■ v„ + ^ + (2) 

where B is magnetic field, and 

y = ap|*p|2 + ^|M/p|4 + a„|vI/„|2 + ^|vI/„|4^c|vI/pnvI/„|2 

(3) 

is the potential term. We begin with a discussion of 
the simplest case of two s-wave condensates (so ~ 
jvE'ple*'*'' and '^n — |^n|e"^" are complex scalar fields 
which discribe proton and neutron condensates corre- 
spondingly). In the above expression 



and 



v„ = (n/2TO„)V(/>„ 



Vp = {h/2mp)V4'p — (2e/TOpc)A 



(4) 



(5) 



are superfluid velocities of neutron and proton conden- 
sates. The key feature of this system is the Andreev- 
Bashkin effect ISj: due to interaction between two su- 
perfluids the particle current of one of the condensates 
is carried by the superfluid velocity of another so the su- 
perfluid mass current of protons and neutrons in such a 
system is liHIT^: 



(6) 



where p^" = p"^ is the superfluid density of one of the 
condensates which is carried by superfluid velocity of an- 
other. Because of the Andreev-Bashkin effect the charged 
supercurrent in this system depends on gradients of neu- 
tron condensate (as it follows from 101, (E))): 
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Let us discuss topological defects, allowed in |(5J), other 
than Abrikosov vortices. 



in the center of the vortex. Let us now impose the fol- 
lowing configuration of (pn- if we go once around the 
vortex core the (/)„ changes 27rn, while if we cover the vor- 
tex loop once in toroidal direction (a closed curve along 
the core) (/>„ changes 2t:1 with n, I being integer. Such 
a situation naturally occurs if a loop is formed around 
rotation-induced vortex line or in case of two interlinked 
loops. This configuration corresponds to a spiral super- 
flow of the neutron Cooper pairs in such a vortex ring. 
Topologically such a vortex is equivalent to knot solitons 
considered in and can also be characterized by a unit 
vector e = (cos (pn sin 9, sin sin 6, cos 9) with a nontriv- 
ial winding. We stress that we do not impose a non- 
trivial winding on (j)p (compare with discussion of knot 
solitons in the two-gap model |lll | where, in contrast, in a 
knot soliton the phases of both condensates must have a 
nontrivial winding number, however, as discussed below, 
neutral-charged mixture is principally different from the 
system in Ref. 0). 

Indeed the nontrivial superflow of neutron Cooper 
pairs induces a drag current of proton Cooper pairs which 
in turn induces a magnetic field which can be calculated 
from 0: 



B ~ curl 



ip cTipP'^mp 
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which can also be written as 
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This self-induced magnetic field gives the following con- 
tribution to the free energy 



^ m, ' ^ ' 
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(10) 



HELICAL NEUTRON VORTEX LOOP 

Let us consider a vortex loop made up of neutron con- 
densate with zero density of neutron Cooper pairs in 
its^core. Let us introduce a new variable 9 as follows: 
— SIT? (^). We will consider a defect where if we 

pPPirin V2/ 

go from the core center to the boundary of the fluxtube 
in a cross section to the vortex, the variable 9 grows from 
to TT. Since at the center of the vortex we have chosen 
that the density of the neutron condensate vanishes then 
indeed there is no drag effect in the center of the flux- 
tube and correspondingly pP" is zero in the core. This 
allows one to chose the boundary condition sin^(0/2) = 



Which is a version of the Faddeev fourth-order deriva- 
tive term analogous to the fourth-order derivative term 
in [ill closely related to the fourth-order derivative 
term in The fourth order derivative terms of this 
type pr ovide stability to finite length topological defects 
UMUM Physically, in a mixture of a charged and a neu- 
tral condensates this effect corresponds to the following 
situation: as mentioned above, the nontrivial configura- 
tion of phase and density of neutron condensate induces a 
charged drag current of proton Cooper pairs which results 
in the configuration of magnetic field This configura- 
tion has the special feature that if the vortex shrinks then 
the magnitude of the self-induced magnetic field grows. 
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We also remark that ffP is a measure of background den- 
sity of proton condensate which is not required to vary 
to produce a knot soliton. 



In the two-gap model in there is a competition of 
the fourth-order derivative term (which corresponds to 
self-induced magnetic field) versus a second-order deriva- 
tive term and a mass term for the third component of the 
0(3)-symmetric order parameter n (the third component 
of n is related to condensate densities in |0 and thus 
it is massive). In contrast, in the present model in the 
competition also participates kinetic energy of superflow 
of neutron Cooper pairs (which is minimized if the vor- 
tex shrinks). A necessary condition for (meta)stability of 
such a vortex loop is that the competition of kinetic en- 
ergy of superflows, gradients of condensate density versus 
the self-induced magnetic field would stabilize the vortex 
loop at a length scale which corresponds to magnitude of 
magnetic field |B(x)| smaller than the field which could 
break proton Cooper pairs. We also emphasis that one 
of the differences with the system of two charged scalar 
fields in is that in the present case the self induced 
magnetic field comes from drag current in the vicinity of 
the core while the superflow of neutron Cooper pairs is 
extended (not localized on length scale shorter or equal 
to p enetration length like the field inducing currents in 
jlllp. We also remark that indeed the effective action 
(0) is assumed being derived from a microscopic theory 
in the approximation of small gradients. Indeed one can 
derive higher-order derivative terms from a microscopic 
theory but this sort of terms, in contrast to the term H10() 
is irrelevant for discussion of the stability of finite-length 
topological defects in this system. Indeed a competition 
between second- and fourth-order derivative terms ob- 
tained in a derivative expansion would stabilize a topo- 
logical defect at a characteristic length scale where all the 
higher-order derivative terms become of the same order 
of magnitude. So, at such length scales the derivative 
expansion fails. We also would like to stress that in the 
present system the knot soliton is prevented against a 
collapse by a finite energy barrier, in contrast to an infi- 
nite energy barrier in the case of the Faddeev's nonlinear 
CT- model considered in mathematical physics That 
is, a zero in proton condensate density, outside core, may 
lead to unwinding of a knot soliton since in a such a 
point the unit vector n is ill-defined and thus the Hopf 
map is ill-defined as well. However the modulus of pro- 
ton condensate order parameter is massive so producing 
such a singularity is energetically expensive. Thus, it is 
a finite energy barrier which prevents a knot soliton in a 
superfluid/superconductor from collapsing. 



AN EXAMPLE OF GENERALIZATION TO 
OTHER PAIRING SYMMETRIES 

Let us generalize the discussion to the case of a mix- 
ture of a spin-triplet neutron condensate and a sin- 
glet proton condensate in order to show that the pic- 
ture does not depend significantly on pairing symmetry. 
The order parameter of the spin-1 neutral condensate is 
|\E'n(x)pCq(x) where {q ~ 1, 0, —1) and C is a normalized 
spinor ' C = 1- Free energy of a neutral spin-1 system 
is (see e.g. (19j): 



Ft = 



2m„ 



(V|*„|)^ 



2m„ 

[co + C2 < S >2] 



*n|'(VC)'-Ai|^r 



(11) 



where < S >= Cq^qj^j is spin. Degenerate spinors are re- 
lated to each other by gauge transformation e"^" and spin 



rotations lA (a, /3, T)=e 



-iFypp-iF, 



', where (a, j3, r) 



are the Euler angles. The topological defects in the neu- 
tral system like this have been intensively studied (see 
e.g. H, ^ charged counterpart of this system in 

ferromagnetic state allows stable knot solitons as it was 
shown in [l^ . 

Let us consider first the ferromagnetic state (which 
emerges when C2 < 0) in context of a mixture of superflu- 
ids. The energy in this case is minimized by < S >^= 1 
and the ground state spinor and density are 19]: C = 
e*('^"-^)(e-*"cos2|, V2cosfsin|,e»"sin2|); |^'„|2 ^ 
— The superfluid velocity in ferromagnetic case is 
n 



v„ = 



2m„ 



[V(0„ — t) — COS pS/a]. So in a mixture of 
a neutral ferromagnetic triplet condensate and a charged 
singlet condensate the equation for charged current is: 
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(^^["^ {<t>n - T)- cos 13V a] 
V pPPrrir, 



(12) 



From this expression we can see that assuming e.g. that 
there is no nontrivial windings in the variables a and /3, 
the system reduces to JTl and thus allows for the topolog- 
ical defects in the form described in the first part of the 
paper. We emphasis that there are no knots of this type 
in a charged ferromagnetic triplet system considered in 
[l^ because in the current equation of a charged triplet 
superconductor, the ratio of the coefficients in front of 
the vector potential term and the gradient term analo- 
gous to V((/()„ — r) does not depend on the carrier density 
and thus one can not obtain a contribution analogous to 
Faddeev term to the free energy by imposing a nontriv- 
ial configuration of the first gradient term in the current 
equation similar to (|12|l in the system jl^]. In a charged 
triplet case the knot soliton may form only as a spin 
texture [l4| . So a neutral-charged mixture with drag 
effect in its magnetic properties is principally different 
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from a genuine charged system. Spin-texture knots can 
be formed in the present system too, as a configuration 
of the order parameter s = (cos (3 sin a, sin /3 sin a, cos (3) 
characterized by a nontrivial Hopf invariant. Such a tex- 
ture generates magnetic field due drag current induced 
by the superflow of the neutron Cooper pairs, which is 
produced by the spin texture. So, in general, there is the 
following nontrivial magnetic energy contribution to the 
free energy functional: 



5127re2 
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(13) 



It must be observed that the spin-texture knot soliton 
is structurally different from the topologically equivalent 
knot of the type considered in the first part of the paper. 
The spin-texture knot is coreless (there are no zeroes of 
the condensate density in it). The third component of 
the order parameter s = (cos /3 sin a, sin /3 sin a, cos /3) is 
massless in this case, thus the spin-texture knot solitons 
in this system are energetically less expensive and have 
larger characteristic size than the topologically equivalent 
knots in the variable e = (cos (pn sin 0, sin 0„ sin 6, cos 6) . 

Let us now consider the "polar" phase of triplet super- 
conductors which is the case when C2 > in The 
energy is minimized then by < S >= 0. The spinor C 
and the condensate density in the ground state are [l9j : 
^ = e'-^" ("Tje-^^sin/?, cos/3, ^e-sin/3); ^ ^/^^^ 

The superfluid velocity in this case is (see e.g. 0|): 
v„ = 2^^*^" which is analogous to singlet case. Thus 
in the antiferromagnetic case the allowed knot solitons 
are equivalent to knots in a mixture of two singlet con- 
densates considered in the first part of the paper. 

We also remark that it is generally assumed that there 
is no proton-neutron pairing in a neutron star becasue of 
large differences in their Fermi energies. 

While we can not make at this stage any definite pred- 
icatoins (which would require large-scale numerical sim- 
ulations), let us however discuss possible mechanisms of 
formation of knot solitons of the discussed above types in 
neutron stars. As it is known, ordinary vortices in super- 
conductors form e.g. as an energetically preferred state 
in external magnetic field. Indeed it is not the only pos- 
sible mechanism of creation of topological defects. For 
instance many defects in liquid helium are created in a 
laboratory without rotation, by various thermal quench 
techniques including neutron irradiation since a symme- 
try breaking phase transition is accompanied by creation 
of topological defects. In case of ordinary Abelian Higgs 
model, created during transition vortex loops shrink, 
while the knot solitons should remain stable. Another 
less common mechanism of formation of vortices is the 
spontaneous vortex state j2oj which emerges when in 



a system coexist superconductivity and magnetism and 
since a vortex carries magnetic field it may have nega- 
tive contribution to free energy functional via Zeeman- 
like coupling terms. So for such systems it is energetically 
preferred to form vortex lattice even without applied ex- 
ternal fields |2oj . Similar mechanism may work in neu- 
tron stars if in a presence of a spin-1 superfluid, there 
is a direct coupling of spins of Cooper pairs to magnetic 
field. In general an effective action of triplet superfluid 
features a Zeeman-like term which is a direct coupling of 
spin S to magnetic field B 



-r]S ■ B 



(14) 



Indeed existence of such a term could result in a range 
of parameters where knots would have a finite negative 
energy if spins of neutron Cooper pairs in the knot soli- 
ton are aligned along the self-induced magnetic field. A 
definite answer to this question is however a complicated 
problem because of competition of many terms. This 
appears being a particularly interesting problem for nu- 
merical simulations. Since in a knot, the magnetic field 
grows if the knot shrinks, it could be that in a such a 
system a formation of a dense ensemble of knot solitons 
is energetically preferred over a spontaneous vortex state 
of Abrikosov vortices. A definite answer to this question 
may however be only obtained in a large scale numeri- 
cal simulation. Thus, if an ensemble of knot solitions is 
formed in a neutron star then one of the apparent conse- 
quences would be its interaction with ordinary columnar 
neutron vortices, then apparently in such a case knot soli- 
tons would disturb a regular lattice of neutron vortices. 



HELICAL PROTONIC FLUXTUBE FORMED 
AROUND A NEUTRON VORTEX 

Above we considered knot solitons which appear due 
to notrivial helical windings of neutron condensate phase. 
In principle there is a theoretical mechamism which 
would allow formation of helical vortex loops of proton 
condensate. Let us show however that a helical protonic 
vortex loop is not a knot soliton and it is not stable. 
Here we stress the most recent studies indicate type-I 
behaviour of proton condensates. Let us now however 
consider the model In that model a neutron star 

possesses a lattice of uniform neutron vortices and a com- 
plicated structure of sparse entangled proton fiux tubes 
(see Fig. 3 in 113 )• the dymanical processes discussed 
in j23| one should expect that entangled proton fiuxtubes 
may dymanically form rings around columnar neutron 
vortices as shown on Fig 1. Let us now consider such a 
ring. The charge current in such a ring is given by eq. l(7|l. 
When we go around such a flux tube, the protonic phase 
(j)p changes 27r, however there is also a current along such 
a vortex due to drag effect by superfluid neutron Cooper 
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A<\>n = 2tz 




A4)p = 271 



FIG. 1: A protonic fluxtube ring around neutron vortex. Due 
to drag by neutron Cooper pairs, the resulting charge current 
in protonic fluxtube is helical. 



pairs which is characterized by a nontrivial phase wind- 
ing of 0„ which changes 2t: when we cover flux tube once 
in toroidal direction. This results in a spiral net charge 
current in such a vortex loop resembling that of a knot 
solition considered in the first part of the paper. Let us 
show however that such a vortex loop is not stable: The 
magnetic field in such a helical fluxtube is given by 



^ ch 
B = —curl 

4e 
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In such a configuration, in spite of helical net charge cur- 
rent, the individual phase configurations of (f)p and </)„ 
are not helical, besides that, the ratio of the vector po- 
tential term to gradient term for i/ip is constant. Thus 
such a helical superflow does not result into a self-induced 
Faddeev-Skyrme-like term, which, if it would be pres ent, 
would significanly affect the considerations in P.l2l|. 



CONCLUSION 

In conclusion we studied topological defects other than 
Abrikosov vortices in an interacting mixture of a neutral 
and a charged condensates. Such a system is believed 
being realized in the interior of neutron stars. We have 
shown that due to Andreev-Bashkin effect the system 
possesses a large variety of knot solitons of different na- 
ture than the knot solitons in the systems studied before. 
We also suggested that due to Zeeman coupling term, 
there could be a theoretical possibility of an exotic inho- 
mogeneous ground state in this system: a spontaneous 
formation of a dense ensemble of knot solitons. 
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